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Abstract. Recent theoretical and experimental progress on the study of ferrotoroidic materials 
is reviewed. The basic held equations are hrst described and then the expressions for magnetic 
toroidal moment and toroidization are derived. Relevant materials and experimental observation 
of magnetic toroidal moment and toroidal domains are summarized next. The thermodynamics 
of such magnetic materials is discussed in detail with examples of ferrorotoidic phase transition 
studied using Landau modelling. Specihcally, an example of application of Landau modelling to 
the study of toroidocaloric effect is also provided. Recent results of polar nanostructures with 
electrical toroidal moment are dually reviewed. 

1. Introduction 

According to the common point of view, ferroelastic, ferroelectric and ferromagnetic materials 
constitute the family of ferroic materials pQ. More recently ferrotoroidic materials have also 
been included in this family Pi- Ferrotoroidics describe materials where toroidal moments 
show cooperative long range order. Ferrotoroidic materials intrinsically belong to the class of 
multiferroic materials |1]. Ferrotoroidicity spontaneously emerges at a phase transition from a 
paratoroidic to a ferrotoroidic phase in which both time and spatial inversion symmetries are 
simultaneously broken ElElilE]. The order parameter for this transition is toroidization. Note 
that in ferroelectrics only the spatial inversion symmetry is broken whereas in ferromagnets 
only the time reversal symmetry is broken. In ferroelastics neither symmetry is broken; only 
the rotational symmetry is broken |5]. 

Ferrotoroidal order can be understood in terms of an ordering of magnetic-vortex like struc¬ 
tures characterized by a toroidal (dipolar) moment. This order is also related to (asymmetric) 
magnetoelectricity, i.e. aij ^ aji [3]. In the present article we are mainly concerned with mag¬ 
netic toroidal moments [21 El [71 El E] as observed for instance, in LiCo(P 04)3 [ID]- Electrical 
toroidal moments can also exist in nanostructures such as polar dots miEiE] but not as a 
long range ordered state in bulk materials as no symmetry is broken. In Section 7 we will sum¬ 
marize recent results on nanoscale ferroelectric materials which exhibit electric toroidization as 
a consequence of dipolar vortex formation. 

Figure [T] shows the symmetry properties of the four ferroic vector order parameters. Po¬ 
larization is a polar vector, magnetization is an axial vector (it contains a sense of time), and 
toroidization is an axio-polar vector. Note that strain is a (second rank polar) tensor order 
parameter. In addition, there are physical properties described by a second rank axial ten¬ 
sor such as magnetogyration HUES] likely present in the spin-half antiferromagnet potassium 
hyperoxide, KO 2 , as well as in CdS, (Gaa;Ini_a;) 2 Se 3 , Pb 5 Ge 30 ii and Bii2Ge02o- 

Therefore, the table can be conceivably generalized to include tensor ferroics. The two 
entries in the left column would be strain (second rank polar tensor) and magnetogyration 


(second rank axial tensor) but at present the tensor analogs of polarization and toroidization 
that would complete the table have not been properly identihed. Note that in principle this 
idea could be generalized further to third (and higher) rank tensor ferroic properties. We intend 
to present these results elsewhere in the near future. 

Any ferroic order is usually accompanied by domain walls [1]. Indeed, ferrotoroidic domain 
walls have been observed in LiCo(P 04)3 using nonlinear optics, i.e. second harmonic generation. 
Litvin has provided a symmetry based classification of such domains [IB] as well as ferrotoroidal 
crystals HZIIIB]. Another material, BCG (Ba 2 CoGe 207 ), also exhibits spontaneous toroidal mo¬ 
ments [T9|. Similarly MnTiOa thin dims show ferrotoroidic ordering PD]. Multiple ferrotoroidic 
phase transitions have been studied in Ni-Br and Ni-I boracites um- Interestingly, some 
quasi-one dimensional materials such as pyroxenes [22] also exhibit ferrotoroidic behaviour. 
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Fig. 1: Symmetry properties of the four vectorial ferroic orders. 


2. Basic field equations 

We will introduce here the toroidic moment and toroidization following ideas published by 
Dubovik et ah in Ref. [S]. Assume distributions of charges p(r) and currents j(r) localized 
in a given region of space. These distributions create electric and magnetic fields that satify 
Maxwell equations. In the presence of matter, within the continuum dipolar approximation, 
these equations are usually expressed as [23] : 


r9D 


= J, 

(1) 

V D 

= P, 

(2) 

(9B 


= 0, 

(3) 

VB 

= 0. 

(4) 


The electric displacement, D, and the magnetic held, H, are dehned as: 
D = £:oE + P, 

H = —B-M, 
ho 


(5) 

( 6 ) 








where P and M are the electric and magnetic polarizations of the medinm, and Eq and /xq 
the electric permittivity and the magnetic permeability of free space. Electric polarization (or 
simply polarization) and magnetic polarization (or magnetization) are introdnced as volume 
densities of the electric and magnetic moments which are dehned from a multipole expansion 
far from charge and current distribution of the electric scalar potential, (p, and the magnetic 
vector potential, A, at dipolar order respectively [2l]. These potentials are dehned from eqs. 
([^ and (|^ as: 

E = -V(p-A, (7) 

B = V X A. (8) 

We will see that the dipolar approximation is not sufficient for some peculiar (electric and 
magnetic) moment conhgurations. In this case, higher order terms in the expansion must be 
taken into account. 

Now let us assume a situation where the conhguration of magnetic or electric moments is 
spiral-like. This is illustrated in Fig. For this conhguration, the magnetization (or polariza¬ 
tion) is along the 2 ;-axis, while the projection on the xj/-plane, which has a circle-like conhgura¬ 
tion, is zero. This circle-like conhguration can be understood as being originated by a toroidal 
conhguration of loop-currents or electric dipoles in the magnetic and electric cases respectively. 
For this kind of conhgurations M or P alone do not provide enough information about the 
ordering of magnetic/electric moments. Irene A. Beardsley [25] noticed that in this case an 
arbitrary amount of a divergence-free magnetization/polarization distribution can be added 
to M or P without ahecting the external held created by the distribution of magnetic/electric 
moments. This divergence-free term can be written as the curl of some vector that characterizes 
the circle-like conhguration of magnetic/electric moments in the xy-plane. In the magnetic case, 
this vector is the magnetic toroidization, V X T^, while in the electric case it is the electric 
toroidization, V x T®. Note that the existence of magnetic toroidization implies that both, 
time-reversal and spatial-inversion symmetries are broken. Therefore, magnetic toroidization 
is represented by an axiopolar (or time-odd polar) vector. In contrast, no broken symmetry is 
associated with electric toroidization. As we will discuss later in Section 7, these toroidizations 
are related to the moment of the distributions of magnetic and electric moments respectively. 
In practice, this means that the magnetization M must be replaced by M -|- V x T^, while in 
the electric case, the electric polarization P must be replaced by P -|- V x T^. 

When electric and magnetic toroidizations are taken into account, the new macroscopic 
Maxwell equations take the same formal expressions as the standard ones after the following 
redehnition of the helds. 


D -)■ D = £oE -h P-h V X T® = £oE-F P, 
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H ^ H = —B - M- V X = —B -M. 
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Therefore, the new macroscopic Maxwell equations read: 
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Fig. 2: Spiral-like configuration of moments along the surface of a cylinder. 


The energy density accounting for the interaction of the generalized polarization and mag¬ 
netization with external helds are given by E ■ P and B ■ M respectively. Hence, the cou¬ 
pling energies of the held with electric and magnetic toroidizations are JE ■ [V x T®]d^r and 
J B ■ [V X T^] d^r respectively. These terms can be written in the form: 


y E • [V X T®](iV = y V ■ (T'^ X E)(iV + /(V x E) ■ 
Jb-IV X T"](iV = ■ (T" X B)(iV + J (V x B) ■ 


(15) 

(16) 


where the hrst terms on the right-hand sides of both eq. (15) and eq. (16) are zero taking into 
account the Gauss theorem. In the electric case, V x E = 0, and thus this energy vanishes. In 
the magnetic case, it is given by J(V x B) ■ T^d^r. Therefore, this shows that the conjugate 
held of the magnetic toroidization is V x B. 

We can now generalize the ideas discussed above and assume systems with spiral-like conhg- 
urations of toroidal moments (either electric or magnetic) arising from toroidal conhgurations 
of electric or magnetic moments. This kind of double-vortex conhguration should be character¬ 
ized by divergence-free vectors expressed as the curl of higher order toroidal moments. These 
moments are denoted as hypertoroidal moments [26]. In the case of magnetism, for instance, 
this means that in the macroscopic Maxwell equations the toroidization T^, should be replaced 
by: 


^ + V X TJ^), (17) 

where is the magnetic hrst-order toroidization and is the magnetic second-order 

toroidization (or hypertoroidization). Hence, magnetization should be replaced by: 

M ^ M(2) = M +V X Tff) + V X (V X t;^)). (18) 

Of course this idea can be formally generalized to any order by dehning higher order hyper¬ 
toroidal moments (see Fig. [^. At order n we will have: 

M(„) = M + V X Tjf) + V X V X Tj^) + ... + V X .... x V x Tjf) 


(19) 






Indeed, one can proceed similarly in the electric case. Therefore, Maxwell equations at the nth 
order, similar to equations PH . can be established by dehning the following helds: 



( 20 ) 

( 21 ) 


With similar arguments as those given above, it is easy to show that the held conjugated 
to the nth order magnetic hypertoroidization is, V x .... x V x B. 
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Fig. 3: Magnetic moment, toroidal moment and the generation of successive hypertoroidal 
moments. 


3. The magnetic toroidal moment and toroidization 

Similar to polarization and magnetization, electric and magnetic toroidizations are dehned in 
the continuum approximation as volume densities of electric and magnetic toroidal moments 
respectively. Since the symmetries associated with electric toroidization are trivial (no change 
of sign is expected either under spatial inversion or under time reversal) no phase transition 
to an electric toroidal phase should be envisaged. Actually, this is consistent with the fact 
that toroidal moment associated with electric moment vorticity is not expected to occur in the 
thermodynamic limit [27]. From the discussion in the preceding section, it seems intuitively 
reasonable to foresee that the toroidal moment should be related to the moment of the dis¬ 
tribution of magnetic moments. This is what we will discuss in the rest of this section where 
we will introduce the toroidal moment based on the multipolar expansion beyond the dipolar 
approximation. We will also introduce a second dehnition based on symmetry considerations 
which is of interest from a more macroscopic thermodynamic point of view. 


In the former case we consider a finite distribution of steady currents J(r). The vector 
potential of this distribution is given by 
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where R is the vector position of a point P, r the vector position of the volume element dv and 
V the volume of the distribution, and the Coulomb gauge (V • A = 0) has been assumed. The 
multipole expansion of A(R) (about r = 0) takes the form [21], 
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It is easy to see that the zeroth-order term vanishes for a steady current distribution for 
which the continuity equation yields V ■ J = 0. The first order term in the expansion is the 
dipolar term that can be expressed as. 
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where m is the magnetic dipolar moment defined as. 
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The next term can be expressed as the sum of magnetic quadrupolar and toroidal contributions. 
The quadrupolar part is given by, 

(A.quad)i^'^ = —^ijkdkNNl^ (26) 

where Eijk is the Levi-Civita symbol and qu is the magnetic quadrupolar moment given by. 
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which is a traceless symmetric tensor. The toroidal term is given by, 

aS = V(t.V)i+«(R), 

where t is the toroidal moment that can be expressed as 
t = -/ rx[rx J(r)](ir’. 


(28) 


(29) 


This pseudovector represents the dual antisymmetric part of the complete tensor which appears 
in the second order term of the multipole expansion. Defining m(r) = x J(r)] as the 
distribution of magnetic moments, the toroidal moment can be written as. 
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( 30 ) 





which indicates that the toroidal moment can be nnderstood as the moment of the distribution 
of magnetic moments (see Fig. 3b). 

For a discrete distribution of N point charges of mass rria localized at positions Fq with 
velocities Uq,, the current density can be written as 

N 

J(r) = - r„). (31) 

a=l 

The magnetic moment then takes the form: 
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where iiIq is the magnetic moment of charge a. Similarly, the toroidal moment of interest here 
can be written as, 


t = ^ ^ ga(r„ X [r„ X Uc,]) = ^ X m„]. 

0=1 0=1 

For a system consisting of a distribution of N spins localized at positions Fq,, iiIq = gfio^a, 
where g is the gyromagnetic ratio and fiB the Bohr magneton. Therefore, the corresponding 
toroidal moment is given by 

1 ^ 

ta = -gHo'^ira X S„) 

0=1 

It is worth noting that treatment similar to the one discussed above can be developed from 
the multipolar expansion of the electric scalar potential (j). In particular, from the first order 
term the electric dipolar moment can be defined as 


(34) 



p = y rp{r)dv. (35) 

y 

Once the toroidal moment is introduced, magnetic toroidization, hereafter simply denoted 
as toroidization, T, is defined as the volume density of toroidal moments. That is T = dt/dv. 
We have already seen that its corresponding conjugated held is V x B and hence, the energy 
density of a distribution of toroidal moments characterized by a toroidization T is given by 
E = -T-{V X B). Therefore, a net toroidization might be induced by means of a current 
density J = V x B. Schmid [2] noticed that reversing toroidal dipoles by means of such a 
held in order to modify toroidization appears unfeasible since this would require the action of 
coherent circular currents of very small size (comparable to the unit cell of the crystal). 

The observation of toroidization in the absence of applied helds indicates the existence of 
long range order associated with toroidal moments. This long range order is usually denoted as 
ferrotoroidic order and should be related to some kind of coupling between toroidal moments. 
Taking into account the basic symmetries of the toroidal moment, the occurrence of ferrotoroidic 
order supposes the simultaneous breaking of spatial inversion and time reversal symmetries. 

Materials with toroidal moments are expected to intrinsically display magnetoelectric cou¬ 
pling. In these systems, an applied magnetic held breaks inversion symmetry and thus induces 
polarization. On its turn, an applied electric held breaks time reversal symmetry and induces 


magnetization. Therefore, we expect that these materials respond to applied electric and mag¬ 
netic helds according to the following equations, 


P = + (36) 

M = a E-t-XmP) (37) 

where Xe Xm respectively, electric and magnetic susceptibility tensors, and is the 
magnetoelectric tensor (all are rank-2 tensors). From a thermodynamic point of view, if the 
free energy of the system is these polarizations and magnetizations should be expressed as 
—dT/dYi and —respectively. Therefore, we expect that the free energy of a magne¬ 
toelectric term is of the type J7n-e = —E B (or —afjEiBj, in coordinate notation). The 
decomposition of this magnetoelectric term into pseudoscalar, vector, and symmetric traceless 
terms enables one to express Em-e in the form (see Ref. m, 

Era-e ~ “E ■ B — T' ■ [E X B] — Qij[EiBj + EjBi], (38) 


where T' is a vector with the same symmetry properties as toroidal moment and toroidization. 
Identifying this vector with toroidization T supposes that its conjugated held is G = ExB. This 
assumption is in agreement with recent experiments that showed that toroiodal moments can 
be controlled by this held [28]. Therefore, polarization, Pt, and magnetization, M^, intrinsically 
associated with the energy term — G-T, induced, respectively, under application of electric and 
magnetic helds, can be expressed as 
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B X T, 
T X E, 


(39) 

(40) 


where we have taken into account that T ■ [E x B] = B ■ [T x E] = E ■ [B x T]. 

In the far-held approximation, from multipolar expansions of electric (scalar) and vector 
potentials, the toroidal held G can be expressed as. 


G = E X B = A{p X m) -|- B{p x r) -|- G(r x m). 


(41) 


where, p and m are electric and magnetic dipole moments, and the coefficients A, B, and C 
are coefficients that decay with distance r as r“®, and r~'^ respectively. Taking into account 


eq. (41), it is worth pointing out that if p and m are parallel, then G = 0 as expected. On the 


other hand, G has maximum strength when p and m are perpendicular. 

Taking G as the held conjugated to toroidization, suggests the following alternative dehni- 
tion of the toroidal moment: 




(42) 


This dehnition neglects residual terms in eq. (41) associated with the magnetic and electric 


moments. The choice is supported by the fact that these terms decay (with distance) faster 
than the magneto-toroidal one. 

It is worth noticing that this dehnition of the toroidal moment is not strictly equivalent to 


the dehnition in eq. (29) resulting from the second order term in the multipole expansion of the 


vector potential. Nevertheless, it is, in fact, expected to provide a good measure of the toroidal 
moment in systems which are simultaneously ferroelectric and ferromagnetic |29|. In these 







systems the coupling of t (or the toroidization obtained as the volume density of this toroidal 
moment) to G leads to a magnetoelectric response similar to that of magnetic ferrotoroidics [3]. 
In general, however, a non-zero toroidal moment should be possible even in antiferroelectric and 
antiferromagnetic systems. Indeed, resonant x-ray diffraction observations of orbital currents 
in CuO provide direct evidence of antiferrotoroidic ordering [30]. Actually, these situations can 
only be considered when the standard dehnition (arising from the multipole expansion) of the 
toroidal moment is taken into account. Note that a multipole expansion including the toroidal 
moment has been considered in |31j . 

4. Materials and relevant experimental results 

At present direct measurements of toroidization or toroidal moment seem very difficult. Present 
experimental techniques can only detect magnetization and magnetic moment, for instance from 
polarized neutron scattering or Lorentz microscopy. In principle, these techniques should be able 
to detect specihc arrangements of magnetic moments characterized by toroidal moments that 
might order to yield net toroidization and thus, ferrotoroidal order. In practice, this appears to 
be unfeasible. 

Indirectly ferrotoroidal order can be inferred from an asymmetric magnetoelectric response. 
Therefore, this needs measurement and analysis of the appropiate magnetoelectric tensor com¬ 
ponents. Sannikov [32] has argued that observation of 7 ^ aji is an indication of possible 
ferrotoroidic order. It is however important to take into account that this condition is not suf- 
hcient to justify the occurrence of toroidization. Asymmetric behaviour of the magnetoelectric 
tensor has been reported for some boracites (G 2 phase of Co-I and Ni-Cl boracites). It has been 
reported also for some oxides such as Ga 2 -xFea ;03 and Gr 203 [331IM] . 

Visualization of ferrotoroidal order requires an experimental technique that is sensitive 
to both space inversion and time reversal broken symmetries which is the inherent feature 
associated with ferrotoroidal order. As shown by Van Aken and co-authors [10] non-linear 
optics offers this possibility. These authors used optical second harmonic generation (SHG) to 
resolve ferrotoroidal domains in LiGoP 04 . Similar experiments were already carried out some 
years before [351136] but were much less conclusive in relation to the existence of ferrotoroidal 
domains. In this technique, electromagnetic light held E(a;) of given frequency is incident on 
a crystal and induces a polarization at double the frequency which acts as a wave source. 
The symmetry affects the corresponding susceptibility. This means that the second harmonic 
generation light from domains with opposite order should have a phase shift of 180°. 

LiGoP 04 crystallizes in the orthorhombic Pnma olivine structure [371 |3H]. If displays unique 
properties including large linear magnetoelectric effect and large Li-ionic conductivity. Go^’*' 
ions belonging to (100) Go-0 layers carry the magnetic moments that are strongly coupled by 
superexchange Go-O-Go interactions. Layers, however, are only weakly coupled by higher order 
interactions. Thus, the system behaves as a magnetic 2-d system to a very good approximation. 

Antiferromagnetic order in the material occurs below a Neel temperature = 21.4 K. Due 
to large magnetocrystalline anisotropy the magnetic moments are conhned to directions lying 
within b-c planes, approximately 4.6° away from the b axis. Actually, the Go magnetic moments 
are not completely compensated, and the system shows a small net magnetic moment which, in 
fact, is not consistent with the orthorhombic symmetry and can only be understood assuming 
a small monoclinic distortion. Interestingly, the monoclinic symmetry allows for a non-zero 
dielectric polarization and a non-zero toroidal moment (along the pseudo-orthorhombic a-axis) 
to occur. 


Van Aken et al. SHG experiments |T0] detected fonr different domain states. It was shown 
later from symmetry considerations [32] that the fonr domains are eqnivalent with different 
orientations of the net magnetic moment. Thns they carry toroidal moments with signs and 
directions mntnally conpled. The fact that the SHG signal intensity is observed to disappear 
precisely at the Neel temperatnre, corroborates the conpling between magnetic and toroidal 
order parameters. 

In addition to Ba 2 GoGe 207 [12] and MnTiOs thin films [201 toroidal moments have been 
considered in BiFeOa and related multiferroics ini im SO]. Another candidate material is the 
magnetoelectric MnPSa as indicated by nentron polarimetry [13|. Based on an initio calcnlations 
the olivine Li 4 MnFeGoNiP 40 i 6 is possibly a ferrotoroidic material H- Note that apart from 
qnasi-one dimensional materials called pyroxenes [22] the toroidal moment in the molecnlar 
context is also of interest, e.g. in dysprosinm triangle based systems [15]. In a related context 
a physical realization of toroidal order is an interacting system of disks with a triangle of spins 
on each disk [16]. 

It is worth noting that the existence of a trne long range ordered ferrotoroidic phase has only 
been established in a snfficiently reliable way in very few cases and, perhaps the clearer evidence 
has been provided by Van Aken et al. SHG resnlts [10]. In other cases, only indirect resnlts 
snggest the existence of snch a phase. The fact that interaction between toroidal moments is very 
weak as indicated by the short range dipolar interaction (see eq. (41)), snggests that any small 
amonnt of disorder in the material is enongh to yield a toroidal glassy state, which represents 
a frozen state with local order only m The existence of toroidal glass in Nio, 4 Mno, 6 Ti 03 
has been foreseen from the behavionr of the magnetoelectric response which was observed to 
strongly depend on cooling history [IH] . This is indeed a very interesting resnlt snggesting that 
materials which are candidates to display ferrotoroidal order shonld also be analysed within 
this point of view. In fact, possible observation of toroidal glass completes the qnartet of ferroic 
glasses, namely spin glass, relaxor ferroelectrics and strain glass 


5. Thermodynamics 


We consider a macroscopic body where the vector ferroic properties, namely polarization, P, 
magnetization, M, and toroidization, T, coexist. Only part of the polarization and magnetiza¬ 
tion will be assnmed to be intrinsic, thns originating from preexisting electric and magnetic mo¬ 
ments. The remaining part will arise from the toroidization originating from magnetic toroidal 
momemts in the presence of external magnetic and electric fields respectively. That is, 

P = P* + Po 

M = M, + Mt 

For this kind of closed systems, the fnndamental thermodynamic eqnation reads 

dU = TdS + ¥.-dP + B- dM, (45) 


(43) 

(44) 


where U is the internal energy density, S the entropy density and r the temperatnre. 


Taking into acconnt eq. (39), the term E ■ dP can be expressed as. 


E-dP = E-dP, + E-dPt = E-dPi + G-dT + T-[Ex dB]. 


(46) 


Similarly, nsing eq. (40), the term B ■ dM can be expressed as. 


B ■ dM = B ■ dMi + B • dM* = B ■ dM^ + G ■ dT + T ■ [dE x B]. 


(47) 





Therefore, 


dU = TdS + E-dPi + B- dMi + G-dT + d{G- T). (48) 

Helmholtz, T, and Gibbs, Q, free energies are dehned as follows, 

jr = U-rS (49) 

g = J^-E-P-B-M = J^-E-Pi-B-Mi-2G-T. (50) 

Their differential expressions are, 

dJ^ = -Sdr + E-dP + B-dM, (51) 

dg = -Sdr -P -dE-M-dB, (52) 

which can be alternatively expressed as, 

djr = -SdT + E-dPi + B-dMi + G-dT + d{G-T), (53) 

dg = -Sdr -Pi-dE-Mi-dB-T -dG. (54) 


Note that this expression suggests that we can assume that the three ferroic properties, po¬ 
larization, magnetization, and toroidization can be assumed as independent (vector) quantities 
thermodynamically conjugated to the electric, magnetic and toroidal helds, respectively. 

The response of the system to applied electric and magnetic helds is given by the generalized 
susceptibility. 
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Diagonal terms, Xe Xmi dehne electric and magnetic susceptibilities. These susceptibilities 
have two contributions. The intrinsic contribution, given by Xg, = dPi/dE (= —d'^g/dE'^), 
and Xm^ — (9Mi/(9B (= —d‘^g/dB‘^), and the toroidal contributions arising from toroidization. 
These last contributions are given by. 
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A toroidal susceptibility can also be dehned as. 


Xt = dT/dG = -d^^g/dG^. 


(58) 


Neglecting the intrinsic contributions to the polarization and magnetization, this toroidal sus¬ 
ceptibility can be expressed as, 
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(59) 


which shows that it is related to Xet ^^d Xmt- 

Cross terms dehne the magnetoelectric coefficients. Maxwell relations require that second 
derivatives of g are independent of the order in which they are performed. Therefore, this yields 

=T = 

a = a. 


(60) 



Assuming that the whole magnetoelectric interplay arises from the existence of toroidization, 
the magnetoelectric coefficient is given by, 

^ dPt 
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-- = 7 — X B + T X I, 
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aM. dExT 


= Ex-+Ix 


(61) 


(62) 


aE aE 

where I is the identity tensor. Notice that thermodynamic stability requires that ^ is positive- 
dehnite. This implies that both Xe Xm must be positive-dehnite and, Xe Xm > « 

Thermal response is determined by the second order derivatives of Q involving temperature. 
On the one hand, the heat capacity C is given by 

d^g C 




(63) 


Taking into account Maxwell relations, the derivatives involving temperature and helds satisfy, 

(64) 


and 
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It can also be obtained that 
d^G d'^G dS 


(9T 


(65) 


( 66 ) 


drdG dGdr 

These expressions determine the cross-response to electric or magnetic held and temperature. 
They are adequate for the study of thermal response of the materials to applied external helds 
which are commonly denoted as caloric ehects. These ehects are quantihed by the entropy 
change that occurs by isothermally applying or removing a given held, and the temperature 
change that results when the same held is applied or removed adiabatically. From a practical 
point of view, materials displaying large caloric ehects are nowadays of great interest thanks 
to their potential use in energy harvesting, and particularly, in refrigeration applications [5U] . 
In general in ferroic and multiferroic materials, large caloric ehects are expected in the vicinity 
of phase transitions to ferroic and multiferroic phases due to the expected strong temperature 
dependence of thermodynamic properties ISH. In the case of toroidal materials, caloric ehects 
have been analyzed from a theoretical perspective in Ref. ^2]. The entropy change induced 
by application of an electric held, (0 —> E), which quantihes the electrocaloric ehect, can be 


obtained from integration of eq. (64) as. 


E 


f dP 

AS(t, 0^E) = --dE. 

J OT 


(67) 


Similarly the entropy change induced by application of a magnetic held, (0 —?■ B), which 


quantihes the magnetocaloric ehect, can be obtained from integration of eq. (65) as. 


A^(r,0 B) = 


dM 

Ih 


■ dB. 


( 68 ) 

















Taking into account eqs. (43) and (44), the preceding entropy changes characterizing elec¬ 
trocaloric and magnetocaloric effects, can be, respectively, decomposed into two terms asso¬ 
ciated with intrinsic contributions and contributions arising from the toroidal moment. The 
intrinsic electro- and magnetocaloric terms are respectively. 


A5i(r,0^E) = 


dr 

(69) 

T‘ • 

OT 

(70) 


A^i(r,0^B) = 


The contributions arising from the toroidal moment can be written in the form. 


ASt{T,0^E) 

ASt{T,0^B) 




dMt f fdT ^ _ 


(71) 

(72) 


A change of entropy can be isothermally induced by application of a toroidal field G = E x B. 
This entropy change characterizes the toroidocaloric effect, and when taken into account with 
eq. (66), it can simply be expressed as, 

G E 


A5(r,0^G) = j^.dG = J 


dT 


dr 


Bx— -dET / —xE -dB 


0 0 
= AS't(r, 0 —> E) -|- AS't(r, 0 —> B), 


9T 


dr 


(73) 


which shows that the toroidocaloric entropy change is simply the sum of the electrocaloric and 
magnetocaloric contributions associated with the toroidal moment as expected. 

Similar expressions can be written for electrically and magnetically induced adiabatic tem¬ 
perature changes. The corresponding total change can be obtained by taking into account that 
from eqs. (63), (64) and (65), the constant entropy condition (adiabaticity in thermodynamic 

(74) 


equilibrium) can be expressed as. 


C , dP dM 

—dr -|- ■ dE -|- —— ■ dB = 0. 

T OT OT 


Therefore, the adiabatic temperature change induced by application of an electric field is given 
as, 

E 

T- r)P 

(75) 


r T dP 

At{S, 0 ^E) = J ^^-dE, 


and the adiabatic temperature change induced by application of a magnetic field as. 


Ar(^,0 B) = 


r aM 
C~^ 


■ dB. 


(76) 














On its turn, the adiabatic temperature change induced by application of a toroidal field is given 

by, 


r dT 


E 


r 


C dr 


C 


dT\ 


T fdT 


At{S, 0^G)= -—■dT= - [Bx — ] ■dE+ / - —xE -dB 


dr 


C \ dr 


(77) 


0 0 0 

where the last two terms in the right-hand side correspond to the sum of electrocaloric and 
magnetocaloric contributions associated with the toroidal moment. That is, 


At{S, 0 ^ G) = AniS, 0 E) + Artis, 0 B) 


(78) 


It is worth noting that both AS{t,0 —)■ G) and At{S,0 —?• G) vanish if either E or B is 
zero or if they are parallel. For the sake of simplicity let us consider that E = [E, 0, 0) and 
B = (0,S,0). In this case, G = (0,0, ES) and assuming electric and magnetic isotropy, 
P = (F, 0, 0), M = (0, M, 0) and T = (0,0,T), and eqs. (73) and (77) are simply expressed as 


EB 


dT 


dT 


B 


dT 


AF(r,0^EF)= / —d(EF) = F / —dE + E / —dB, 


dr 


dr 


dr 


(79) 


and 


EB 


B 


Ar(S,0^ EB) = / If = B IE ^eJ if dB. 

0 0 0 

respectively. 


(80) 


6. Landau and Ginzburg-Landau modelling and domains 

Phase transitons in ferroic and multiferroic materials are associated, as we have discussed 
above, with some symmetry change. This change is captured by an order parameter which 
is zero at temperatures above the transition and non-zero below it. Landau and Ginzburg- 
Landau theories provide reliable expressions of the free energy of the materials in the region 
of the transition in homogeneous and non-homogeneous cases, respectively. The approach is 
phenomenological in nature and its combination with the thermodynamics formalism provides 
a powerful method to study macroscopic and mesoscopic behaviour of ferroic and multiferroic 
materials. This approach enables one to relate measurable quantities to the input parameters of 
the theories that can be determined either from experiments or from first-principle calculations. 
In Landau theory the free energy is expressed as a series expansion of the order parameter. As 
this free energy must be invariant under the symmetry operations of the system, only those 
terms allowed by symmetry are included in the series expansion. 

In ferrotoroidic materials toroidization is the primary order parameter but magnetization 
and polarization must also be included in the free energy. So far, few Landau models have 
been proposed to account for ferrotoroidal transition in specihc materials. Sannikov [2l] al¬ 
ready proposed a model to account for the anomalous behaviour of the component 032 of the 
magnetoelectric tensor near the cubic (43ml') to the orthorrombic (m'm2') phase transition in 
boracites. 

Based on a group theoretic analysis Sannikov has provided a free energy for ferrotoroidic 
phase transitions in boracites [2III32]. It consists of a usual double well F(T) = aT‘^ + bT‘^ in the 




toroidization T and harmonic terms cP^ and in polarization and magnetization. In addi¬ 
tion, it has symmetry allowed conpling terms between varions components of the three dipolar 
vectors, namely of the form PiTj, MiTj and the trilinear conpling PiMjTk. The free energy is 
then analysed to obtain a phase diagram in terms of the free energy coefficients indicating the 
varions ferrotoroidic phase transitions. Note that Litvin has extended the symmetry analysis 
for ferroics to ferrotoroidic materials including the possible toroidic domains and domain walls 
[isiiniiiH]. This analysis is very helpful in obtaining the free energy for ferrotoroidic phase 
transitions for any crystal symmetry. 

Ederer and Spaldin, taking into account that the symmetries which allow for a macroscopic 
toroidal moment are the same that give rise to an antisymmetric component of the linear 
magnetoelectric effect tensor, proposed the simplest possible Landau free energy that describes 
a phase transition between a paratoroidic and a ferrotoroidic phase that includes the energies 
associated with the effect of electric and magnetic held on polarization and magnetization, and 
their coupling to toroidization [53]. It has the following form, 

F(r,T.P.M) = - r“)r2 + tcr‘+ 

- B-M-E-P + kT-(P X M), (81) 


where Xp and Xm are electric and magnetic susceptibilities respectively, Aq is the toroidic 
stiffness, and C > 0 is the nonlinear toroidic coefficient, k measures the strength of the mag¬ 
netoelectric coupling. The last term in the previous free energy represents the lowest possible 
order coupling term between the three order parameters consistent with the required space and 
time reversal symmetries. 

Minimization of the free energy (81) with respect to polarization and magnetization provides 
the equilibrium values of polarization and magnetization. They are given by the following 
equations. 


P = XpE - X T), (82) 

and 

M = - Xmfi;(T X P). (83) 


We can now solve these two equations assuming (for simplicity) that E = (E,0,0) and B = 
(0,5, 0) and, assuming isotropy, P = (P, 0,0), M = (0,M, 0) and T = (0,0, T). We obtain: 

P = XpE - KXpXmBT + 0[T‘^) ~ XpE - aB (84) 

and 

M = XmH - KXpXmET + O(r^) ~ XmB - aE, (85) 


where we have neglected the nonlinear magnetoelectric effects in the above two equations. The 
magnetoelectric coefficient a = nXpXmT is a quadrilinear product of electric susceptibility, 
magnetic susceptibility, the coupling constant k and the toroidization. Thus, either for k = 0 
or r = 0 there is no magnetoelectric effect. 


Substitution of P (84) and M (85) in the free energy (81) gives the following general type 


of effective free energy: 

3 = FoiE,H) + ^Ao{T-T,)T^ + ^(3T^ + ^CT^ + XT, 


P 


( 86 ) 






where, 


Fo = 

-liXpE^ + XmB^), 

(87) 

Tc = 

^2 

Tc + -^XpXmiXmB^ + XpE^], 

(88) 

/S = 

^i^^XmXlEB = Sn^XmXpK 

(89) 

A = 

f^XmXpEB. 

(90) 


The effective free energy (86) corresponds to the free energy of a toroidal system subjected to 
an effective applied toroidal field A (proportional to G). Interestingly, the coefficient of the third 
order term, /3, depends also on A (and thus, on G). When G = 0, (3 = 0 then (86) describes a 
paratoroidal-to-ferrotoroidal second-order phase transition. Under the application of a toroidal 
field G ^ 0 (and (3 ^ 0), the transition becomes first-order for |A| > |/3pC'^/27. It is worth 
pointing out that the addition of nonlinear terms in eqs. (84) and (pSj) would lead to higher 


order terms in the expansion (86) that go beyond the minimal model. However, within the 
spirit of the Landau Theory, it is expected that such terms are not essential. 



Fig. 4: Toroidal order parameter as a function of temperature (in arbitrary units) for three 
different values of the coupling parameter k = 0.90, k = 1 and k = 1.05 and selected values of 
the applied toroidal field G. The arrow in the lower panel indicates the occurrence of the first 
order transition. Taken from |52] . 


It is worth noting that the free energy (86) does not include a term directly coupling G 
and T, which is in agreement with the thermodynamic formulation developed in section ^ In 
Ref. [53], however, P, M and T were assumed to be independent order parameters and a term 







































-G T was included in the free energy. The same effective free energy rtSdh is obtained, but in 


this case the parameter A is given by, {i^XmXp ~ 1)EB. Similar results are obtained with this 
model which leads to a richer physics due to the fact that the A and (3 terms can have different 
sign. The model has been used in [52] in order to study toroidocaloric effects in ferrotoroidic 
materials. The entropy of the system can be obtained as. 


S(r,T,G) = 


dF 1 


(91) 


where T{t,G) is the equilibrium value of the toroidal order parameter which is a solution of 
dFe/dT = 0. Then, the change of entropy isothermally induced by application of a toroidal 
held is obtained as. 


S{T, G = EH) - S{T, G = 0) = --Ao[T\T, G = EH) 
It is easy to show 


T\T,G = 0)]. 


(92) 


l^see [ni]) that this expression coincides with the general thermodynamic 
expression (73). The variation of toroidization as a function of temperature for three different 


values of the coupling parameter (k) is shown in Figure 4. The corresponding isothermal entropy 
change (AS) or the toroidocaloric effect is depicted in Figure 5. 



Fig. 5: Toroidocaloric effect, i.e. the isothermal change in entropy, as a function of transition 
temperature (in arbitrary units) and selected values of the applied toroidal held G for three 
diherent values of the coupling parameter k = 0.90, k = 1 and k = 1.05. Taken from |52j . 

The present Landau approach to ferrotoroidic materials can be generalized by including 
symmetry allowed gradient term (VT)^ [9], i.e. the Ginzburg term. This should allow one to 
study domains and domain walls in ferrotoroidic materials as observed in LiCo(P 04)3 using 


















optical second harmonic generation techniqnes [TO]. With doping induced disorder in such 
materials we expect that novel phases such as toroidic tweed and toroidic glass should also 
exist and remain to be observed experimentally with certainty jlH]. With symmetry allowed 
coupling of strain to toroidization, if we apply stress to such a crystal we expect toroidoelastic 
effects, i.e. a change in toroidization with hydrostatic pressure or shear. We expect that these 
important topics will be explored in near future. 

7. Electric ferrotoroidics at the nanoscale 

We have already discussed in Section 3| that no broken symmetry is associated with electric 
toroidization which is consistent with the fact the formation of electric moment vortex is forbid¬ 
den in the thermodynamic limit. However, the situation can drastically change when the scale 
of the material decreases towards the nanoscale. It has been predicted that vortex structures 
can be stabilized below a certain critical temperature in both ferroelectric and ferromagnetic 
nanodots [SllESI. These zero-dimensional structures have been studied in detail from ab initio 
simulations of an effective hamiltonian under appropriate boundary conditions m- They nu¬ 
merically simulated Pb(Zr,Ti )03 ferroelectric nanodisks and nanorods under open circuit-like 
electric boundary conditions and found vortex states comprising electric dipoles that form a 
closed structure yielding spontaneous electric toroidal moment. Indeed, in recent years electro- 
toroidic behavior has been found experimentally jSSl EH ESI ESI ED] • Atomistic simulations of 
KTaOa seem to indicate an incipient ferrotoroidic response wherein quantum vibrations sup¬ 
press the formation of polar vortices [6T]. In addition, electrogyration or electric held induced 
optical rotation has been predicted in materials exhibiting electrotoroidic behavior [I3]. They 
suggest that these nanostructures are potentially interesting for data storage applications. Note 
that the effect of long-range elastic interactions on the electric toroidization in a ferroelectric 
nanoparticle has been considered in [62] . 

From a practical point of view, using these nanostructures based on electric toroidal moment 
as writing memory nanodevices is not straightforward since, in principle, these moments cannot 
be switched by standard methods as they are unaffected by applied electric helds. However, 
several solutions have been envisaged. In [63] it has been shown that vortices, both electric 
and magnetic, can be manipulated by inhomogeneous static helds. The coupling with elasticity 
enables controlling vortices by mechanical load, which gives rise to a rich temperature-stress 
phase diagram [M]. Another possibility, considers nanodots with the shape of nanorings with an 
oh-central hole. The interest in these objects relies on the fact that they are characterized by a 
transverse hypertoroidal moment, which is a polar vector and thus sensitive to an homogeneous 
applied electric held [65] . Following similar ideas, it has been recently proposed [66] and numer¬ 
ically predicted that ferroelectric nanotori can possess an homogeneous hypertoroidal moment 
as well as exhibit the coexistence of axial toroidal moment and hypertoroidal moment phases. 
In these nanoscale objects the hypertoroidal moment could be manipulated by an homogeneous 
applied electric held. Note that as a diherent application, metamaterials based on the electric 
toroidal moment have been proposed |67j . 

8. Conclusions 

With the recent emergence of magnetoelectric and multiferoic materials the fourth primary 
ferroic property, namely ferrotoroidicity, has gained special attention. In the present article we 
have developed a general thermodynamic framework for the study of phase transitions, domain 


walls and caloric effects within the context of Landan theory in ferrotoroidic materials snch as 
LiCo(P 04 ) 3 . Both magnetic and electric toroidal moments were considered although the latter 
can only exist in polar nanostructures. The generalization to hypertoroidal moments was also 
presented. We discussed a variety of materials where toroidic order has been observed. In the 
presence of sufficient disorder the other three primary ferroics exhibit glassy behaviour, namely 
as spin glass, relaxor ferroelectrics and strain glass [19]. Similarly, toroidal glass has been poten¬ 
tially observed as well in the study of dynamics of a linear magnetoelectric Nio, 4 Mno. 6 Ti 03 |15] . 
Presence of toroidal moments is also an indicator of magnetoelectric coupling in the material 
|68j . Similarly, toroidal magnon excitations in multiferroics relate to magneto-optical effects 
[69]. In this article we did not consider the helds and radiation from moving toroidal dipoles 
which is also an important area of research [TQlE]. Clearly, exploration of toroidal phenomena 
is a fertile area of research with a great potential for both fundamental science and device 
applications. Indeed, toroidal metamaterials have been proposed, studied and experimentally 
realized [72l [73l [71| including in double-ring [75] and double-disk m structures. 
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